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^ 1. Introduction 

q The ultradiscrete KdV equation 

2 Till + T t J = max (T* +2 + T l ]+l - 1, T* +1 + T/+ 1 ), (1) 

is obtained from the bilinear form of the discrete KdV equation [I] by a limiting 
procedure called "ultradiscretization" [2]. By setting 

U* = Tj +1 + Tj +1 -T}+l-Tj, (2) 

this equation is transformed into 

^ +1 = min(l-C/*,£(^-^, +1 )), (3) 

j'=jo 

where jo is determined by the boundary conditions. This equation is known as the time 
evolution rule of the Box and Ball System (BBS) [3], which is a cellular automaton with 
soliton like behavior in spite of the simple time evolution rules. The boundary conditions 
which are actively studied are Uj = for \j\ ^> 1 (infinite systems) and Uj +L = Uj for 
some L G Z >0 (periodic systems). We may choose j as jo — — 00 for infinite systems 
and can also set jo for periodic systems under some conditions. 

These systems have good mathematical structures as well as continuous and 
discrete ones. For infinite systems, in previous papers jl], [S], we proposed a recursive 
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representation which corresponds to the notion of vertex operators. As an analogue of 
determinant-type solutions, the ultradiscretization of signature-free determinants (called 
"Permanent s" ) is discussed in [6] and the relationship between this type of solution 
and ultradiscrete soliton equations is discussed in [7] and [8]. Approaches revealing 
combinatorial properties of the solutions are presented in [H] , [TO] and [TT] , by expressing 
them as maximum (minimum) weight flows of a planar graph. 

For periodic systems, algebro-geometrical methods are considered to be most 
suitable and many topics are studied in the sense of this method. The initial value 
problem of periodic BBS (pBBS) is solved bypassing the analysis for "discrete" (non- 
ultradiscrete) elliptic curve in [T2] . An ultradiscrete closed method however is presented 
in [13]. A direct correspondence of these methods is presented in |14j . 

A description of the dynamics of the BBS using the representation theory is 
presented in [15]. This approach is applicable to both infinite systems and periodic 
ones. 

In this paper, we first consider a discrete quadratic function with a parameter and 
discuss the properties of this function, especially the values of the dependent variables 
where it attains its maximum. We prove that this maximum is a solution of the 
ultradiscrete KdV equation. We propose some examples of such solutions, which include 
well-known pseudo-periodic solutions and soliton solutions of the equation. Finally, we 
discuss the recursive representation we proposed earlier. The approach used here does 
not depend on algebro-geometric methods, but a few results of discrete convex analysis 
are used. 



2. Discrete quadratic form 



Let N be a natural number, T>i (i = 1,...,N) be a discrete interval [a,i, b{] C Z, a 
discrete semi-infinite inteval (— oo,6j], [ai, oo) or Z and T> — T>\ x . . . x T>^ C Z* N . We 
consider a discrete quadratic function for m e D with parameters z e WL N , defined as 

U ' ' (4) 



/(z; m) = -VrMm + ^zm, 
where the matrix A G mat(M, N) is given by 



(A) 



i-1 



-L i + 2^Vt l + 2{N-i)Qi 



i.k 



1=1 



-2Qi 

and parameters L, 6 R satisfying the relations: 

i < n a < n 2 < . . . < n N 

i-1 

2^Qi + 2{N-i + 1)0^ < Li (i = 



(< = *) 

(z < k) 
(i > k) 



N). 



(5) 



(6) 
(7) 
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By introducing a matrix M 6 mat(R, N) given by 



(8) 
(9) 
(10) 



and n e 7L N expressed as 

m = M^t, 
the function / is transformed into: 

0(z;A*) = ^Bfi + ^Mfj,, 
where B = f MAM, is found to be 

' (A) M 
(^4)i-i,i-i + (A) M + 4fi i _ 1 



It should be noted that these two quadratic functions are equivalent because M -1 is 
expressed as 

' 1 (i > k) 
(i<k) 

Proposition 1 B is negative definite. 



i — k — i) 

i = k>2) 
i — k — 1) 
i - k = -1) 
\i-k\> 2) 



(m- 1 ; 



(12) 



Proof We separate the matrix B into B = B\ + B2, with Si and B2 expressed as 



(Bi)i, 



Lk 



i = k = 1) 
i = k>2) 
i-k — 1) 
i-k = -l) 
\i-k\> 2) 



(A) ltl +2Sl 1 

(A)^!,^! + 2Qi_i + (A) M + 2ft; 

— (-^.Ji-l.i-l — 2ftj_i 
2ftj 



B 2 = 2diag(-fti, Six — ft 2 , . . . , ftjv-i - ftjv). 
Here, 5i can be transformed into 

XM-^B.M- 1 = diag((A) 1>1 + 2ft x , (A) N>N + 2Q N ), 
which is negative by condition (|7fi and B 2 is at least non-positive by (j6|. 

Due to this proposition, the quadratic function /(z; m) = g(z; /x) always has a 
maximum and the set on which it elements attains its maximum is finite. 



(13) 

(14) 

(15) 
□ 



Fact 2 The quadratic form ^/iBfi satisfies the relation: 



where A and V means 

(pAp')i := min(/ii,/i-) 



(/x V //); := maxQLti,/4). 



(16) 
(17) 
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This fact is a famous result of convex analysis and a proof is presented, for example, 
in EE61. 



Definition 3 We introduce an ordering for fi by 
H > yl > fjj. forl<i<N. 



Proposition 4 The set of fi which realizes the maximum of (10) has a unique maximum 
element with respect to this ordering. 

Proof We assume the existence of two different local-maximal elements fi and fi' and 
denote /j," = n A y! and y!" = y V //. Due to Fact [2] and identity: y + y! = y" + n'", 
we obtain 

^yBy+'zMn+^n'Bfi'+'zMfi' < B /j," '+ i zMij!'+^ n By n, +kMiJi' n . (19) 

However, this inequality is actually an equality because y and y' yield a maximum, i.e., 
\x" and ft'" also yield a maximum. This contradicts that y and /j,' are local-maximum 
because y!" > n, fi'. □ 

Definition 5 For each t,j G Z ; Zei 

2*=ffi-jl + C (20) 
n = i (fi 1 ,...,fi JV ) (21) 
C = t (C 1 ,...,C JV ) (22) 

1 ='(!,- ••,!)• (23) 
VFe define [A * as ^ which yields the maximum of g{z t ,\y) and which is the maximal 



element for the ordering (18) for such fi. We also denote ni^ := Mfij'* and T^' 1 
as the maximum value of g(zp fi) = /(z*;m) , i.e., 

if* = max /(zj; m) = /(*}; mf >•*) = ^zj; M f >•') (24) 

J m&Z N J J J J J 

We simply denote Tj and m* when we do not need to consider the value of N . 

Let rriN be the iV-th component of m and m be a subvector consisting of the first 
N — 1 components of m, / (z*-; m) is transformed into 

N-l 

/(z*; m) = - ( - L w + 2 ^ fi,)m^ + z N m N + /(zf 2 ^; m), (25) 

i=i 

where / is the discrete quadratic function for m written in 

/(z; rh) = ^VhArh + 'zm, (26) 

matrix A is a submatrix of A consisting of the first N — 1 rows and columns and 
~t-2m N ^ jjN- 1 j s a subvector of z*~ 2mjv consisting of the first N — 1 rows. In other 
words, A is obtained by replacing L, — >■ Lj — 2fi; (i = 1, . . . , N — 1) in the definition 
p\ for N — 1. It should be noted that the condition (frh is also satisfied for / when all 
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parameters Lj and Qi are fixed. Therefore, / has the maximum value. By denoting this 



maximum for z* as T- 1 , we obtain the recursive form 



= max \\ ( - L N + 2 V + W + ff -^"'j . (27) 

Because of the above discussion, the first N — 1 components of are equal to 

m g Z^ -1 which yields the maximum of f(z t ~ mN ; m), where is the iV-th component 
of We also define / = for iV = 1 (i.e. T^' 1 = 0), for consistency. 

3. Behaviour of maximizing vectors 

To prove that T^'* solves the ultradiscrete KdV equation, we study the behaviour of 

(N),t 

m) . 

Lemma 6 /t*- satisfies the relation: 

H) < < /tf 1 . (28) 

Proof We denote z* = z, /x* = /x, fx t j _ 1 = n', /i" = /i A fx' and fj,'" = /i V /i'. Due to 

Fact [2] and z*_ x = z* + 1, we obtain 

-fyBn + *zMli + -VS^i' + \z + l)Mp' 

< ^li"Bli" + bM/x" + )^p!"Bn'" + *(z + l)M/t"' - lM(/i'" - ft'). (29) 
Here, one has H.M(fx'" — fx') > because of the definition of ft'". Thus, we obtain 

^V*£/i + ^Mli + -W + *(z + l)Mfjf 

< ^li"Bli" + feMji" + ^WB/j!" + \z + l)M/t"'. (30) 

By the same discussion as in the proof of Proposition |4| li" and li'" must yield a 
maximum for and (t,j — 1) respectively. By virtue of the maximality of li' and 
the definition of li'", we obtain ft!-_i — fx' — lx'" > /x = /t*-. 

The proof for ft* < /t^ is completely the same. □ 

Theorem 7 m* -2 the same as m* or expressed as m* = m*~ 2 + ej 6 P /or some 
1 < i < N . Here, ej is t/ie i-t/i canonical basis vector. 

Before giving the proof, we prepare a lemma. 
Lemma 8 Under Theorem^ one has 

T (N),t + 2 + T (N),t-2 < ^rp(N),t + 2Qn _ (31) 
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Proof By denoting '* +2 = m, m^'* -2 = m' and z* = z, one has 

T ( N ),t+2 + T (N),t-2 = It m ^ m + t( z + 2 ri) rn+ IWylm , + '(z-20)m , .(32) 

Due to Theorem [7j m — m' is equal to 0, or + e^. In the case of m = m', m* is 
also equal to m. Then, we obtain 

r N)>t+2 ( Ar)>t _ 2 = 1 h mAm + t zm 

3 3 2 2 

= 2T\ NU < 2T\ n) * + 2Q N . (33) 

In the case of m = m' + e^, we obtain 

rp(N),t+2 rp(N),t-2 
3 3 

= -*(m' + e.j)A(m / + e*) + *z(m' + e*) + - m'Am' + *zm' + 2f^ 
< 2Tf ] ' 1 + 2fii < 2T J (JV) '* + 2Q N . (34) 
In the case of m = m' + ej + e^, due to identity: 

-\m + ei)A{m + e,) + -\m + e k )A(m + e fc ) 

= + e, + e fe )A(m + e; + e fe ) + |Wm - (A) itfcj (35) 



we obtain 



rp{N),t+2 + = l t(m , + e . M(m , + e?) + tz(m , + e . } 

+ ^ < (m / + e k )A(m + e fc ) + *z(m' + e fe ) + 2fi, + 2fi fe + (A\ k . (36) 



Here, by the definition of (A)i k , one has 



i-1 



2a+2n,+U)„: { tT 1 (37) 

2^max(i,fc) 7^ ^0 

and this value is less than 2Qn, in both cases, by virtue of ([6| and (J7|). □ 
Now, let us prove Theorem [7j 

Proof We employ the inductive method for N. It is clear that (m*~ 2 )jv < (m*) jy 
because ix l ~ 2 < /a*-. By denoting (m*-)jv = m N and (m* -2 )^ = m! N and employing (27), 
Tj '* and T- '* 2 are expressed as 

AT-l 

= \ ( - L N + 2 £ a) m% + z N m N + ff~^^ (38) 
i=i 

jV-1 

rp(N),t-2 ^^- Ln + 2J2^<+ (zn ~ 2n N )m' N + jf (39) 



i=l 
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and by the maximality for T^- '* and Ty^'* -2 , one h a s 

N-l 



T J (N) ' t >U-L N + 2j2Vi)(m N -l 

8=1 



+^(m iV -l)+ff- 1) ''- 2m - +2 (40) 



j 

N-l 



T^ t - 2 > 1 -(-L N + 2j2^)(m' N + l) 2 

i=i 

+ [z N - 2Sl N ) (m' N + 1) + f\ N ~ l) ^-\ (41) 



Adding these equalities and inequalities, we obtain 

N-l 



iv — i 

> (l n -2^0^ {m N - m' N - 1) - 2fi 7 
i=i 



_l_j,(AT-l),t-2mjv+2 _|_ y(A r -l),t-2m^ v -4 _ j,(N-l),t-2m N _ rp(N-l),t-2m' N -2 

By assuming that this theorem holds for N — 1, we can apply Lemma [8] repeatedly and 
obtain 

rp(N~l),t-2m N +2 Sp(N-l),t-2m' N -4 _ rji(N-l),t-2m N _ rp(N-l),t-2m' N -2 
3 3 3 3 

>-2Q N (m N -m' N -2). (43) 



Therefore, (42) is transformed into 

N 



> (L N -2^VL^j(m N -m' N -l). (44) 

i=l 

For N = 1, we also obtain the same result because Tj ^'* is equal to 0. To satisfy this 
inequality, itln — m' N must be less than 1 because of ^ and Q. 

Thus, we obtain (m*-)jv - 1 < (m* -2 )^ < (m$) N . If (m*-)jv = (m'" 2 )^ + 1, by 



the relation (27), the i-th components of m* and m* 2 are the same for i < N — 1. If 



(m*-)jv = ( m j~ )jvj it is clear by the assumption of the induction. □ 
Lemma 9 m* +1 is equal to m* or express as m* = m*- +1 + e.j 6 "D /or some 1 < « < AT. 
Proof Note that ^* = ^*~ 2 or 

i=i 

for some i by virtue of Theorem [7j Due to Proposition |6j can be equal to 

expressed as 

i 
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for some io > i or 

ii hl-i i 

i'=i 2 +i i'=i 2l +i i'=i 

for some io > i\ > z 2 > . . . > %n-\ > %n > i- 



We have to prove A*j+i cannot take the form (47). We assume A*$+i takes the form: 

i\ i 

A*j+i = A*5- X] e i'"X ei '' (48) 
i'=i 2 +i i'=i 

which is the simplest case of ( |47j ) and equivalent to m*- +1 = — e i2 + — e^. Let 

n 

a*' = a*5 - Yl ei ' ( m ' = m 5 ~ e ^ + e n) ( 49 ) 



i'=ta+l 



and 



a*" = a4 - X ei ' ( m " = m 5 - ei ') • ( 5 °) 

i'=l 

We also denote A^+i = A*'" ( m j+i = m '") and /x*. = /j, (m*- = m). It should be noted 
that (A*)i a = (a 4 ')^ = (a 4 ")^ = (a 4 "')^- We let mi and be the vectors consisting 
of the components z 2 to iV of m and m'". We also let m 2 and m 2 be the vectors 
consisting of the components 1 to 22 — 1 of m and m'". In other words, m = < ( < mi, *m 2 ) 
and m'" = *(W 1 ,W 2 ). Then, m' and m" are also expressed as m' = *(W 1 , i m 2 ) and 
m" = *(*mi, Vn 2 ). 

By employing (25) repeatedly, one can rewrite 

/(z*; m) = S(m 1} t) - fj^j + Jtf-^ ; m 2 ) (51) 

/(**; m') = 5K, t) - + 7(z^ i2 ; m 2 ). (52) 

Here, S is a function which depends on mx 6 7j N ~ 12 and t but not on j, and / depends 
on z e W 2 and m 2 G Z 22 . The maximality of m leads to 

S(mi,t) - S^m^t) > 0. (53) 

However, by the same discussion for j + 1, one has 

/(z$ +1 ; m") = S( mi ,t) - ^ 2 (j + 1) + /(z^ 2 ; m' 2 ) (54) 

/(z$ +1 ;m"') = 5(mi,t) -^(7 + 1) + 7(<i 12 ; m' 2 ). (55) 

Since m" attains the maximum value of /(z* +1 ; ■ ■ ■) and the maximum element for the 
ordering (18), we obtain 

S(mi;t) -Sim'^t) < 0, (56) 

which is a contradiction. It is clear that we can extend this proof to the general case of 
{471 □ 
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By employing all these lemmas and enumerating all possibilities, we finally obtain 
the important theorem: 



Theorem 10 Vectors m* +2 , m*^ 2 , m j +i > m j+i' m j an ^ m j+i are expressed as 



,t+2 



t+1 



(mf 2 , m;+ 2 , m* +1 , m*+i, m*, m$ +1 ) 



m* +1 , m* +1 , m* +1 , m* +1 , m* +1 , m$. +1 ) + A 



Lj + l, ±±± j + l 



i+n 



where A zs one of the following: 

(0,0,0,0,0,0) 
(e;, 0,0, 0,0,0) 
(ei,ei,d,d,0,0) 
(ei,ei,d + e fc ,d,0,0) 
(ei,ej,d + e fe ,d + e fe ,0,0) 
(ej, ej, ej, ej, ej, 0) 
(2ej, ej, e^ ej, ej, 0) 

(^i &h ^i; ^i; ^i; ^i; 0) 

(ej + e fc , ej, d + e k >, d + e k >, e k , 0). 
i/ere ; 2 and k satisfy i > k and d is expressed as 



e io + . . . + e 



*2Z-1 



/or some i > i\ > i<i > . . . > 121-1 > in > k' > k > 0. 



4. The ultradiscrete KdV equation 



Theorem 11 The function T satisfies the ultradiscrete KdV equation: 



rpt+1 _j_ rpt 



-3 
max 



(T' +2 + T* 1 -1,t; +1 +T;+ 1 ) 



(57) 

(58) 
(59) 
(60) 
(61) 
(62) 
(63) 
(64) 
(65) 
(66) 

(67) 



(6£ 



Proof We denote m* +1 = m and z* = z. It is sufficient to prove each case of A in 
Theorem [lOl 

(I) In the case of A = (0, 0, 0, 0, 0, 0) 



rpt-\-2 I rj-if 



tmAm + %z + 2f2 - l)m + -kirn 



?m 



- Wlm + *(z + f2 - l)m + -Wm + '(z + f2)m 

2 2 



T 



t+i 



(69) 



1 



7j+ 2 + 71 = -Wm + *(z + 2fi)m + -Wm + *(z - l)m 

J J ) ) 



rpt-\-2 j_ Tit . rpt+2 

1 j ^ J i 



rpt 



(70) 
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(II) In the case of A = (e i; 0, 0, 0, 0, 0) 

T t+2 + T t > v m + ei )A(m + e*) + *(z + 2fi - l)(m + e*) + -Wm + k 
J 2 2 

= -*(m + ej)A(m + e;) + %z + 2f2)(m + e*) + ^kAm + \z - l)m - 1 
= Tj +2 + Tj +1 -l. (71) 

it+2 _|_ rpt rpt+\ _|_ rpt- 



The proof for T*+ 2 + Tj = T*\ + T m is the same as that for (K39h by the definition 



of Tj+i- 

(III) In the case of A = (e^, e i; d, d, 0, 0) 

T'Xl + T\ = Am + ei)A(m + e 4 ) + *(z + 2fi - l)(m + e,) + -Wm + fern 
J 2 2 

1 1 

= -*(m + ei)A(m + e») + *(z + 2f2)(m + e*) + -Wm + *(z - l)m - 1 

2 2 

= T'+ 2 + T* - 1 (72) 



T;+ 2 + Tj > + d)A(m + d) + *(z + 2fi - l)(m + d) 

+ ^'(m + d)A(m + d) + fe(m + d) 

= -\m + d)A(m + d) + *(z + fl - l)(m + d) 
2 

+ J'(m + d)A(m + d) + '(z + n)(m + d) 

= Tj+i + Tj +1 (73) 
(IV) In the case of A = (e i; e*, d + e^, d, 0, 0) 



TjXl + Tj> -%m + d + e k )A(m + d + e k ) + \z + fi)(m + d + e fc ) 



+ + d)A(m + d) + \z + n - l)(m + d) + (fi* - 1) 

> + r; +1 (74) 

The proof for Tjg + Tj = Tj +2 + Tj +1 - 1 is the same as g. 

(V) In the case of A = (e*, e,, d + e^, d + e^, 0, 0) 
The proof is completely the same as for (III). 

(VI) In the case of A = (e i; e i; e», e,, e*, 0). 
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T^Xl + Tj = -\m + ei)A(m + e*) + \z + 2ft) (m + e,) 
+ ^\m + e,)A(m + e 4 ) + \z - l)(m + «*) 

> ^*( m + ei)A(m + ei) + *(z + 2fi)(m + e*) + ikk + *(z - l)m 



= Tj+ 2 + T/ +1 > + T| +1 - 1 
The proof for T|+ 2 + Tj = Tj^l + Tj +1 is the same as g. 
(VII) In the case of A = (2e i; e,, e^, e,, e^, 0). 



(75) 



The proof for Tg( + Tj = T^i + T] +i is the same as (69) and by virtue of (35), 
one has 

" J 2 

I rvi 



7?g + = + e,)A(m + e<) + *(z + 2ft - l)(m + «*) 



+ -*(m + ei)A(m + e*) + *z(m + e; 



1, 



m + 2e, ; )A(m + 2ej) + (z + 2f2)(m + 2e 



+ -Wlm + *(z - l)m - 1 + (A) ifi - 2Vt t . 
Due to relation ([7]), we obtain 

TlXl + Tj > T*+ 2 + Tj +1 - 1. 
(VIII) In the case of A = (e, + e^, e», e$, e*, e,, 0). 



(76) 
(77) 



The proof for T*+ 2 + Tj = Tjg + T* +1 is the same as (69) and by virtue of (35), 
one has 

Tj+l + Tj = -\m + ei)A(m + e*) + \z + 2Q - l)(m + e,) 
+ ^( m + e,)A(m + e 4 ) + <z(m + e 4 ) 



1, 



m + 2ei)A(m + 2e t ) + \z + 2ft) (m + 2e0 



m 



(z — l)m — 1 + {A)ii — 2f2j. 



Due to relation ([7]), we obtain 

rpt+2 _|_ T-I* ^ T-it+2 I T-it _ 1 

(IX) In the case of A = (e$ + e*., ej, d + e^, d + e^, e^, 0) 



(78) 
(79) 



By virtue of (35) and i > k, one has (A)^ = — 2f2& = —2Klek- Then, 



it t 



T*%i + Tj = -\m + ei)A(m + <*) + %z + 2ft - l)(m + e,) 
+ ^'( m + e fc )A(m + e fc ) + *z(m + e fc ) 
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= -%m + e { + e k )A(m + e t + e k ) 

+\z + 2f2)(m + e» + e fe ) + ^mAm + \z - l)m - 1 
= Tj +2 + Tj +l -l. (80) 
We also obtain T*ff + Tj > Tj+{ + Tj +1 by the same proof in the latter part of (III). □ 

5. Examples 

In this section, we study several aspects of the behaviour of the functions we proposed. 
We employ the dependent variable Uj defined in (|2) in the plots, because this variable 
is best suited for observing the behaviour we are interested in. In the following figures 
we depict the j- lattice by a row of boxes, containing a ball when Uj is equal to 1, and 
empty whenever Uj = 0. 

5.1. Infinite Domain 

In this subsection, we treat the case where D = Z* N and assume that all parameters (flj, 
Ci, Li) are integers. Before arguing the general case, we first consider a special case 
where the solutions are reduced to well-known ones. Let L\ — . . . = = L, and we 
obtain the identity: 

f(z] +L ; m - 1) = /(«*; m) - \z] - (81) 

because of Al = —LI. Due to this identity, if m yields the maximum of /(z; ■), m — 1 
also yields that of f{z — 1; •). In particular, by substituting z = z*, one has: 

m) +L = m* - 1 (82) 

- 1 

t Ul = T h Y,( m > - ^ + c *) + 2 NL - (83) 

1=1 



By virtue of (83), we obtain the relationship Uj +L = Uj, where Tj expresses a state of 
the pBBS. Indeed, this is known as a standard form of solutions for the pBBS. We note 
that the condition ^ simplifies to 2 YliLi ^« < A which is a famous requirement in the 
analysis of pBBS. Figure [T] depicts an example of such solutions. 

We now consider a generalizations of (81 )-(83 ). If we can find K > and L\, . . . , L^ 



satisfying 

An = -Kl (84) 

for given n e Z^ , we obtain the relationship 

m l j+K = m* - n (85) 
Uj +K = Uj, (86) 

by employing the similar argument. In this case, Tj also expresses a state of the pBBS. 
The main difference from the standard form is that the each block of balls parametrized 
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Figure 1. Uj for N = 2, V = Z 2 , < t < 4, < j < 34, fl x = 1, f2 2 
17, Li = 12, L 2 = 12. The area between dashed lines is one period. 



2, Ci = 22, C 2 



□ 



Figure 2. By plotting C7j for iV = 2, 2? = Z 2 , < t < 4, < j < 34,^ = 1,0 2 
2, Ci = 22, C 2 = 17, Li = 8, L 2 = 12. The area between dashed lines is one period. 



with Ci) emerges (n)j times in a single period. For such cases, we should take N 
to be the number of apparent blocks when we employ the standard form. However, 
in our representation, we may employ smaller N. Therefore, the function is in fact 
a "compressed" representation for this system. Figure [2] depicts an example of such 
solutions. 

Finally we consider the case where L\ , . . . ,Ln are given, as shown in Figure [3] 
We observe that each block of balls parametrized by Ci) has its own "pattern", 
depending on Lj. To discover its global behaviour, we consider the equations: 

(An) 1 = ... = (An) N . (87) 

Since all coefficients of A are integers, by the assumption, we can obtain n e 1, N and 
K > satisfying the relationship: An = —Kl. Therefore, by virtue of the same 
arguments as above, solutions are always periodic, for general parameters. In fact, the 
solution depicted in Figure [3] has a period of 38. We note that it is hard to predict 
the period, directly from parameters. All solutions depicted in Figures [T] [3] take the 
same parameters (fl, Ci) and L2. However, we observe that just a little change for L\ 
provokes a large difference for the periods. 
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□ 



Figure 3. By plotting Uj for N = 2, V = Z 2 , < i < 4, < j < 34,fii = l,fi 2 
2, Ci = 22, C 2 = 17, Li = 10, L 2 = 12. 



Figure 4. By plotting C/J for TV = 2, 2? = [-1, 2] X [-1, 1], < t < 4, < j < 34, fi t 
1, « 2 = 2, d = 22, C 2 = 17, Li = 8, Za = 12. 



5.2. Finite Domain 

We next consider the case where all Pj are finite. Figure |4j depicts an example of such a 
case. Blocks of balls emerge several times but disappear in distant sites. Therefore, this 
solution satisfies the boundary condition Uj = for |j| ^> 1 and solves (|3j) for j = — oo, 
which is the time evolution of the standard BBS. Especially in the case D$ = [0, 1] (for 
all i), each block appears only one time and the solutions express the well known soliton 
solutions for this system. This type of solution is also a compressed representation for 
the state with regularly-positioned blocks of balls. However, this compressed solution is 
easily rewritten from the non-compressed (well-known) one, i.e., it is not a new solution. 

Of course, we can consider the case where the Di are semi-infinite. For such 
cases, the corresponding block appears in a pattern for sufficiently large j but never 
for sufficiently small j (or vice versa). We can also take D to be a product of finite 
interval and infinite one, which expresses the state where some blocks appear infinitely 
many times but other blocks appear only finitely many times. 

Finally, we note that limiting the domain corresponds to truncating summations 
for discrete systems and we stress that the truncated solutions of discrete integrable 
equations cease being solutions. 
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6. Concluding Remarks 

In this paper, we have discussed properties for a class of multi-variable quadratic 
functions and we have proven that these functions solve the ultradiscrete KdV equation. 
We have also proposed a new type of solutions by restricting parameters, for example, 
multi-periodic solutions. 

In our previous papers jl], [5], we proposed a recursive representation of soliton 
solutions for the ultradiscrete soliton equations including the ultradiscrete KdV 
equation. This representation can be considered as a transformation from a soliton 
solution to another one, that is, a transformation between two states in the same 
dynamics — the standard BBS, with infinite and open boundary condition. However, 



the recursive representation (27) cannot be regarded as transformation between two 



"standard" forms of solutions to pBBS, because of the shift of the parameter Lj between 
rp(N),t j,(n i),^ j n Q.jjjgj. worc l S; when we obtain a standard form by applying the 



recursive form (27) repeatedly from the vacuum solution tJ ^'* = 0, the intermediate 



states are not standard forms except for some special cases (for example when all f2j are 



equal). It should be noted that relation (27) is an extension of a recursive representation 



of soliton solutions. When restricting Dn = [0, 1], we can omit the term of m 2 N in (27) 
by replacing Cn — > Cat + L N — 2 ^f^ 1 f2j because of m N = m 2 N for m N = 0, 1. The 
representation presented here is suitable for only analyzing the standard BBS. 



As seen in the proof of Theorem 11| it is important to consider the state of m*, 



which corresponds to each state of the BBS. Especially, the case (IX) corresponds to the 
interaction of solitons. In [11] , it was sufficient to consider only two cases, interacting or 
not, to describe the dynamics of the ultradiscrete Toda molecule equation. However, for 
the ultradiscrete KdV equation, we also have to consider an additional case — "injecting" 
balls. The reason for this complexity results from the introduction of coordinates j for 
boxes to represent the Box and Ball dynamics. 

It is also known that the BBS has waves called "backgrounds" which can take 
various value and travel at speed 1 and arbitrary initial states consist of solitons and 
backgrounds. Cauchy problems can be exactly solved by virtue of an ultradiscrete 
analogue of the inverse scattering method [17] . However, our approach in this paper 
cannot be applied to states that include backgrounds, as it depends strongly on good 
combinatorial properties of the solitons. 

We also note that the discussion in the previous sections becomes much easier in 
the case Di = [0, 1] (for all i) because we consider only limited cases. Specifically, we 
do not need to employ Proposition [TJ or most of Theorem [7j Therefore, we can present 
an induction free discussion, such as proposed in [TTj . 



We believe that discrete convexity plays a very important role for the r-functions 
of the ultradiscrete systems. In this paper however, we only used convexity for some 
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propositions which are not so essential. The quadratic function ^ is certainly convex 
and has good combinatorial properties. However, we cannot obtain good properties of 
the solution Tj f itself. We also note that the independent variables t and j in (24) can 
be extended to real values, but m in Q can not. The discreteness for m is considered 
to be essential. It is expected that a direct relationship between these notions and 
fundamental properties of integrable systems such as the Plucker relations, and can be 
expressed in the language of discrete convex analysis. 
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